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Abstract. The first, self-consistent calculations of the cosmological H2 dissociating UV background produced during the 
epoch of reionization (EOR) by the sources of reionization are presented. Large-scale radiative transfer simulations of 
reionization trace the impact of all the ionizing starlight on the IGM from all the sources in our simulation volume down 
to dwarf galaxies of mass ~ IO^Mq, identified by very high-resolution N-body simulations, including the self-regulating 
effect of IGM photoheating on dwarf galaxy formation. The UV continuum emitted below 13.6 eV by each source is then 
transferred through the same IGM, attenuated by atomic H Lyman series resonance lines, to predict the evolution of the 
inhomogeneous background in the Lyman- Werner band of H2 between 11 and 13.6 eV. 
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SUPPRESSION OF MINIHALO POP III 

STAR FORMATION BY H2 

DISSOCIATING UV BACKGROUND 

Simulations suggest that first stars formed inside mini- 
halos of mass M ^ lO^M© at z ^ 20, when H2 cooled 
the primordial, metal-free halo gas and gravitational col- 
lapse ensued. Formation of first stars inside minihalos, 
however, is suppressed when H2 cooling is suppressed 
due to a strong H2 Lyman- Werner(LW) band radiation 
field. This process occurs when the LW intensity /lw ex- 
ceeds the threshold LW intensity (/lw) threshold, which is 
usually set by requiring that cooling time equals the dy- 
namical time (e.g. Haiman, Rees & Loeb 1997). 

It is important, therefore, to correctly calculate the rise 
of H2 dissociating background in the Universe. Earlier 
estimates (e.g. Haiman, Abel & Rees 2000) found that 
sources of reionization made /lw > (-^lw) threshold long 
before reionization is complete, thus sterilizing miniha- 
los before they could contribute significant reionization. 

The LW background intensity would show spatial 
fluctuation due to the inhomogeneous distribution of 
reionization sources. However, previous calculations are 
based upon a homogeneous universe approximation. In 
these calculations, the sources and the IGM opacity 
were both uniformly distributed, in space with uniform 
emissivity given either by analytical approximation (e.g. 
Haiman, Abel & Rees 2000), or by sum over sources 
found in small-box simulations, too small to account for 



large-scale clustering of sources or follow global reion- 
ization (e.g. Ricotti, Gnedin & Shull 2002; Yoshida et al. 
2003). 

Here we present the first self-consistent radiative 
transfer calculations of the inhomogeneous LW back- 
ground produced by the same sources which reionized 
the universe in a large-scale radiative transfer simulation 
of reionization. We show that the fluctuation is signifi- 
cant enough to induce widely-different evolutionary his- 
tory in the first star formation in different parts of the 
Universe. 



A METHOD TO CALCULATE THE 

INHOMOGENEOUS LW BACKGROUND 

DURING EPOCH OF REIONIZATION 

Calculating the inhomogeneous build-up of the LW 
background is computationally challenging. First, the 
mean free path for LW photons (~ 100 cMpc) is much 
larger than the mean free path for H ionizing photons. 
We need to account for sources distributed over large 
volume and look-back time. Second, LW band photons 
redshift and get attenuated by H atom Lyman series res- 
onance lines as they travel through the expanding uni- 
verse. This requires a multi-frequency radiative transfer 
calculation in a cosmological volume (^ (lOOcMpc)^), 
which is very expensive computationally. 

We have developed a simple, yet accurate method 






. r 1 1 1 1 1 1 1 1 1 1 1 1 r 1 1 1 1 1 1 1 1 1 1 1 J 






3 1 




























■o 0.8 




























^ 0.6 














1 0.4 


L 










J 
















2 0.2 


















, . . 1 1 . > 


, , 1 


, 




1 - 



11 11.5 12 12.5 13 13.5 

E (eV) 




20 40 60 ao 

[(l + z)/(l+20)]°- X r^ 



FIGURE 1. "picket-fence" modulation. A source observed at different places suffers different amounts of attenuation, dilution, 
and redshift. As continuum (flat- spectrum assumed) photons redshift, all the photons which have redshifted to H atom resonance 
lines are turned into low-energy photons. This leaves gaps in observed spectrum, which depend on the distance to the source. The 
modulation factor (right panel) is the effective transmission rate, which drops as the distance increases, until source is completely 
attenuated at rcMpc = 92 [(1 +z)/21]~^-^. With this modulation factor, we are able to compute LW intensity without costly multi- 
frequency calculation. 



to calculate the radiative transfer of LW band photons, 
which alleviates the need for a multi-frequency calcula- 
tion. This is achieved by using a pre-computed modula- 
tion factor, expressed in terms of the comoving distance 
between a source and an observer, that accounts for the 
attenuation of LW band photons from a single source. 
We call this a "picket- fence" modulation factor /mod, af- 
ter the shape of the LW intensity Jy which is composed 
of both attenuated and unattenuated parts (Fig. 1). The 
picket-fence modulation factor is computed by averag- 
ing a normalized Jy (at rcMpc(^s, ^obs)) over frequency in 
the range of [n.5 - 13.6] eV. 

At a given point in space, one then looks back in time 
to search for sources located inside the LW band horizon 
(~ 100 cMpc). For each source found this way, a "bare" 
flux is obtained from its luminosity Ly and the luminosity 
distance, which corresponds to a flux in an optically thin 
limit. This is then multiplied by /mod to account for the 
attenuation by H atom Lyman series resonance lines. 
Since /mod is a function of the comoving distance only, 
an expensive multi-frequency calculation is not required. 
When looking back in time for contributing sources, 
we draw past light-cones in order to account for the 
retarded time effect due to finite light-crossing time (Fig. 
2). Sources are distributed in a discrete sense in time, 
following our reionization simulations results (e.g. Iliev 
et al. 2007). 



INHOMOGENEOUS LW BACKGROUND 

In order to calculate inhomogeneous LW background, 
we use "self-regulated" reionization simulation results 
by Iliev et al. (2007). Their radiative transfer calcula- 
tion is based upon an N-body simulation (1624^ particles 
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FIGURE 2. Conformal spacetime diagram, to be used to cal- 
culate source contribution from cosmological distances. Each 
source contributes LW flux of f^Q^L/D^, where L is luminos- 
ity, Dl is the luminosity distance to the source, and /mod is the 
picketfence modulation factor which accounts for IGM opacity. 
This scheme naturally accounts for the finite light-crossing time 
of radiation. 



and 3248^ cells) result which provides the source cata- 
logue for H-ionizing radiative transfer calculation. This 
resolves halos of mass M/Mq ^10^, thus accounting 
for all atomic cooling halos in the (50cMpc)^ volume. 
Source formation is self-regulated: formation of sources 
inside small-mass halos (10^ < M/Mq < 10^) is sup- 
pressed if H II regions overtake their formation sites. 
We then use a coarse-grained mesh (203^ cells) to cal- 
culate the fluctuating LW radiation field. To account for 
the large LW band horizon, we distribute sources period- 
ically around the box where the LW background is cal- 
culated. 






FIGURE 3. Global history of cosmic reionization and LW intensity evolution. (Left): The accumulated number of ionizing 
photons per baryon vs. the number of dissociating photons per baron. (Middle): Global evolution of ionization fraction < x>. 
(Right): Global evolution of LW intensity J21. Note that if reionization source properties change, especially in ionizing photon 
escape fraction (/esc) and ratio of dissociating/ionizing photons, J21 may differ even at the same global ionized fraction. 



We chose an illustrative case, WMAP3 f2000_250S, 
among the reionization simulation results by Iliev et al. 
(2007). For small mass halos (10^ < M/Mq < 10^), Pop 
III high efficiency emitters with top-heavy IMF are as- 
sumed, with fy = f^' /esc 'Ni = 2000, and with Ni/N^iss — 
15, where /* is the star formation efficiency, /esc is the 
ionizing photon escape fraction, Ni is the number of ion- 
izing photons emitted per baryon, and A/diss is the number 
of H2 dissociating photons emitted per baryon. Forma- 
tion of sources inside large-mass halos (M > lO^M© ) 
are not suppressed even inside H II regions. Pop II low 
efficiency emitters with Salpeter IMF are assumed, with 
fr = /* -/esc -Ni = 250, and withM/A^diss ^ 1. 

Fig. 3 summarizes our result on the global ionized 
fraction < x> and the mean LW intensity < /21 >• We 
find that < /21 > is dominated by small mass halo con- 
tribution initially, while at later epoch by large-mass halo 
contribution in this case. By the time large mass halos be- 
come important, most of small mass halos are already in- 
side H II regions and get sterilized. If reionization source 
properties were different from our illustrative case, how- 
ever, especially in /esc and M /A/diss, the same < x > 
would not mean the same /21 . in general. 

A significant inhomogeneity in the LW background is 
apparent, with fluctuation scale of a few - ~ lOcMpc 
(Fig. 4). This is due to a highly-clustered source distri- 
bution, which shows a fluctuation scale of a similar kind. 
This will induce an inhomogeneous LW feedback, hence 
inhomogeneous minihalo star formation rate. The most 
pristine environment will be found where the most ac- 
tive minihalo star formation activity occurs. Those sites 
would exist where the LW intensity becomes minimum, 
as shown in Fig. 4. These first- star forming regions dur- 
ing the epoch of reionization might be detected by up- 
coming telescopes such as the James Webb Space Tele- 
scope (JWST). 




FIGURE 4. Patchy reionization and patchy H2 dissociating 
background in (50cMpc)^ box at z = 15.7. Contours of thick 
colored lines represent different /21 contours (orange - J2i = 
0.1; blue - J2i= 1), and the black contours represent the ion- 
ization fronts. Minihalos are subject to spatially- varying LW 
feedback effect, thus yielding spatially-varying minihalo star 
formation rate. JWST may able to detect clustered minihalo 
population at z ~ 15 in regions of least LW feedback. 
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